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ABSTRACT 

VERITAS has been monitoring the very-high-energy (VHE; > 100 GeV) gamma- 
ray activity of the radio galaxy M87 since 2007. During 2008, flaring activity on a 
timescale of a few days was observed with a peak flux of (0.70 ± 0.16) x 10~ n cm~ 2 s -1 
at energies above 350 GeV. In 2010 April, VERITAS detected a flare from M87 with 
peak flux of (2.71 ± 0.68) x 10~ n cm~ 2 s _1 for E > 350 GeV. The source was observed 
for six consecutive nights during the flare, resulting in a total of 21 hr of good quality 
data. The most rapid flux variation occurred on the trailing edge of the flare with an 
exponential flux decay time of 0.901q 22 , days. The shortest detected exponential rise 
time is three times as long, at 2.87lg]gg days. The quality of the data sample is such 
that spectral analysis can be performed for three periods: rising flux, peak flux, and 
falling flux. The spectra obtained are consistent with power-law forms. The spectral 
index at the peak of the flare is equal to 2.19 ± 0.07. There is some indication that 
the spectrum is softer in the falling phase of the flare than the peak phase, with a 
confidence level corresponding to 3.6 standard deviations. We discuss the implications 
of these results for the acceleration and cooling rates of VHE electrons in M 87 and the 
constraints they provide on the physical size of the emitting region. 

Subject headings: galaxies: individual (M87, VER J1230+123) - gamma rays: galaxies 



1. Introduction 



M 87 is a giant radio galaxy located in the Virgo cluster at a distance of 16.7 Mpc (IMei et al 
20071). It is believed to harbor a supermassive black hole of mass (3.2±0.9)xl0 9 M^ (jMacchetto et al 
19971 ). derived from gas kinematics, or (6.6 ± 0.4) x 1O 9 M (|Gebhardt et allboilh T derived from 



stellar kinematics. Its jet is misaligned with the line of sight; this, along with the prox imity of 
M 87, allows for detailed observations of its structure in the radio (e.g.. ICheung et al.ll2007r ). optical 



(e.g.. lBiretta et al.lll999l ), and X-ray (e.g.. iMarshall et all 12002 : IWilson fc Yandl2002r) wavebands 



Appa rent superluminal motion is observed in the radio and optical wavebands (jBiretta et al.lll995l . 
19991 ) . Month-scale flaring activity has been obs erved in various energy ranges at the nu cleus and 
at HST-1, the jet feature closest to the nucleus (jPerlman et al.ll2003l ; lHarris et al.ll2009i ). The jet 
knot HST-1 is located 0.85 arcsec (« 69 pc projected) from the nucleus and is resolved from the 
nucleus in the radio, optical, and X-ray energy bands. 

Very-high-energy (VHE) gamma-r ay emission from M 87 was first detected by HEGRA in 
1998/1999 at energies above 730 GeV (jAharonian et al .1 120031 ) and has since been confirmed by 
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H.E.S.S. (jAharonian et alJl200fl ). VERITAS (jAcciari et alJl2008h . and MAGIC ([Albert et aUEOOa ). 
The first gamma-ray flaring activity from M87 was reported by H.E.S.S. in 2005, with the flux 
varying on a timescale of days. The angular resolution of current ground-based gamma-ray in- 
struments is not sufficient to distinguish the different morphological features in M87, which is 
therefore detected as a point-like source in VHE gamma rays. However, given the short timescale 
of the flare, the characteristic size of the gamma-ray emitting region (or of moving regions of low 
gamma-ray opacity between us and the emitting region) is constrained by the light crossing time of 
these features and the relativistic Doppler factor of their motion in the observer's reference frame. 
Under the size constraint, the two most likely region s for gamma-ray pro duction are the unresolved 



nucleus and the HST-1 knot ( Aharonian et al. 



2006 



Cheung et al.ll2007h . 



During the 2005 gamma-ray flare ob served by H.E.S.S., Chandra reported historically maximal 
flaring from HST-1 (Ha rris et al.ll2006l ). Through the causality argument, the timescale of the 
enhanced TeV emission implies an emission region size of ab o ut R < 5 x 10 15 5 cm, where S is 
the Doppler factor of the radiating region. lAharonian et al.1 (|2006l ) preferred the nucleus over 
HST-1 as the VHE gamma-ray production region due to an unrealistically small opening angle 
(~ 1.5 x 10~ 3 5 deg) required to channel energy from the central object to the HST-1 knot. However, 
the Very Long Baseline Array (VLB A) imaged compac t knot s in HST-1 that are not resolved 
with semi- minor axes < 5 x 10 17 cm (jCheung et al.l 120071 ). and lStawarz et al.l (j2006l ) proposed jet 
reconfinement at the HST-1 location which can in turn produce TeV emission. Therefore, HST-1 
remains a candidate for TeV emission. However, a VHE gamma-ray flare in 2008 coincided with 
the historical maximal X-ray flux from the nucleus detected by Chandra, while HST-1 remained in 
a low state at that time and its X-ray flux was below that of the nucleus. In addition, increasing 
radio flux from the nucleus, but not fr om the jet, was obs erved by the VLBA, lasting up to two 
months past the VHE gamma-ray flare (jAcciari et al.ll2009l ). The 2008 observations therefore favor 
the nucleus as the origin of the VHE gamma-ray emission. 

After the launch of the Fermi Gamma-ray Space Telescope in the summer of 200 8, M87 was 



also d etected in the MeV-GeV energy range by the Fermi Large Area Telescope (LAT; lAbdo et al 
20091 ). However, no significant flaring activity was detected in 2009 at any wavelength. 



M 87 has been monitored every year in VHE gamma rays since 2003 by at least one of the three 
major atmospheric-Cherenkov telescope arrays — H.E.S.S., MAGIC, and VERITAS. In 2010, VHE 
flaring activity up to 20 % of the Crab Nebula flux was detected from M 87 in the span of several days 
(jOng &; Mariottill2010l ). and gamma-ray, X-ray, optical, and radio observations were subsequently 
triggered. Detailed results from the VERITAS observations are presen ted in this paper, and th e 



multiwavelength light curve will be presented in a separate publication (jAbramowski et al.1 12011 
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2. Observations and Analysis 



VERITAS is an array of four 12 m diameter imaging atmospheric Cherenkov telescopes located 
at the Fred Lawrence Whipple Observatory in southern Arizona, 1.3 km above sea level. The 
telescopes are situated approximately 100 m apart, forming a convex quadrilateral. Each telescope 
is equipped with a camera of 499 photomultiplier tubes (PMTs) arranged in a hexagonal lattice 
covering a field of view with a diameter of 3.5°. The array is sensitive to photons with energy 
from ~ 150 GeV to more than 30 TeV, with an angular resolution of ~ 0.1° and an effective area of 
~ 10 5 m 2 at 1 TeV. Fur ther description of the V ERITAS observatory and its performance are given 
in lPerkins e^all <|2009h and bolder et al.l (|200(tI ). 



M87 was observed between 2009 December and 2010 May for 53.1 hr. Observations were 
conducted at a range of zenith angles between 19° and 40°, with low elevation excursions (up to 
60° from zenith) during the nights of April 9 through 11 when episodes of flaring were detected. 
More than 95% of the data were taken with the full four-telescope array and the remainder with 
a three-telescope sub-array. To enable simultaneous estimation of source and background signals, 
the data were accumulated in "wobble mode" for which the source is offset from the camera center 
by 0.5° in alternating directions every 20 minutes. The analysis presented in this paper is based 
on 44.6 hr of live time which satisfied data quality and integrity selection criteria. 



The data are analyzed with the algorithm described in lAcciari et al.1 (120101 ). Atmospheric 
gamma-ray shower images are first corrected for dispersion in PMT gain and timing using infor- 
mation obtained from nightly laser calibrations (jHannal 120081 ) . Then, an image-cleaning process 
is applied to select pixels with a signa l significantl y above the night-sky background level. After 
cleaning, the images are pa rameterize d (Hi 



the stereoscopic technique (jHofmann et al 



^ 1985 ) and the shower direction is reconstructed using 



[1993). if at 

least three camera images pass selection criteria optimized for a source with 1% of the Crab Nebula 
flux. The results reported i n this paper have all been confirmed by an independent secondary 



analysis package described in iDaniell (|2008l ) 



3. Results 

During the six-month observation period, M87 was detected at a level of 25.6 standard de- 
viation (a) above the background, with an average flux of (5.44 ± 0.30) x 10~ 12 photon cm -2 s _1 
at energies above 350 GeV, equivalent to 5% of the Crab Nebula flux above 350 GeV. The follow- 
ing sub-sections first present the daily light curve obtained over six months of observation, then 
the April flaring episode light curve binned in 20 minute intervals, followed by the timescale and 
spectral analyses of the April flare. 
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3.1. Daily Flux Over a Six- month Period 



Figure Q] shows the daily flux recorded by VERITAS between 2009 December and 2010 May. 
Applying a constant-flux fit to the daily light curve gives a x 2 /dof value of 269.4/29, a strong 
indication the flux was not constant during the observation period. 

In 2010 February, the MAGIC Collaboration re ported an incre ased activity of M87 with more 
than 10% of the Crab Nebula flux on February 9 (|Mariottil [2OI0I ). At that time, VERITAS ob- 
servations were hampered by poor weather conditions, but M87 was detected by VERITAS in 
a typical state two nights after the MAGIC alert. In 2010 April, VERITAS detected M 87 with 
an elevated flux during a wee k of observations betw een April 5 and 11 and triggered subsequent 
multiwavelength observations (jOng &: Mariottill2010l ). 



3.2. The 2010 April Flare 

Figure [2] shows the light curve binned in 20 minute intervals during the flare. Observations 
around the peak of the flare were carried out up to high zenith angles; this reduces the sensitivity 
at low energies, and as a result, the data points taken at the end of April 10 and the beginning of 
April 11 have larger uncertainties. The flaring episode began with increasing flux during the nights 
of April 5 and 6, reaching 10% of the Crab Nebula flux on April 8. On the following three nights, 
M87 was observed for more than five hours each night. The average flux on April 9 and 10 was 
15% of the Crab Nebula flux, reaching as much as 20% of the Crab Nebula flux in individual 20 
minute bins. The average flux on April 11th was 5% of the Crab Nebula flux. VERITAS continued 
to monitor M87 for two hours each night from April 12 to 15, when the flux level returned to a 
few percent of the Crab Nebula flux, comparable to the low-state flux measured in the past. All 
flux comparison with the Crab Nebula is at energies above 350 GeV. 



3.2.1. Flux Variability Timescale Analysis 

Using the data from April 9 and 10 (MJD 55295 and 55296) when maximal activity occurred, 
we searched for variability within each night. On April 9, fifteen 20 minute exposures were taken 
in total and a constant-flux fit yields a x 2 /dof value of 9.3/14 and a corresponding x 2 probability 
of 0.81. On April 10th, twenty-one 20- minute exposures were taken and the constant-flux fit gives 
a x 2 /dof value of 19.8/20 and a corresponding \ 2 probability of 0.47. In order to investigate 
variability within a single day of observation in more detail, the wavelet analysis described by 



Price et ajj (|201ll ) is applied to the April 9 and 10 data sets. The highest confidence level for the 
April 9 data set is obtained for a variability timescale of 80 minutes. However, the confidence level 
is only 86.2%, implying that on an 80 minute timescale, the evidence for variability is only at the 
level of 1.5 a. The highest confidence level for the April 10 data set is obtained for a variability 
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Fig. 1. — Daily light curve of M 87 observed by VERITAS in 2010. Clear evidence of flaring activity 
is seen in 2010 April (MJD 55291-55298). Trigger alerts sent by MAGIC on February 10 (MJD 
55237) and by VERITAS and MAGIC on April 9 (MJD 55295) are indicated by vertical lines. 
The average nightly flux during the peak of the flare exceeds 10% of the Crab Nebula flux at the 
same energy threshold of 350 GeV. A constant spectral index of 2.5 was assumed for the daily flux 
calculation. 
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Fig. 2. — VERITAS light curve with 20 minute binning during the flare period between 2010 April 
5 and April 12 (MJD 55291-55298). The flux scale is the same for all six panels, and dashed lines 
indicating 5%, 10%, and 20% of the Crab Nebula flux are included. A constant spectral index of 
2.5 was assumed for the flux calculation. 
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timescale of 160 minutes at 97.5%, or 2.2 a. Therefore, no evidence for intra-night variability is 
found. 

Figure [3] shows the daily light curve of the April flaring episode. To characterize the timescales 
of the flare, an exponential function of the form <3? = po e ( i_55290 )/pi j s fitted to different periods 
of the April daily light curve by x 2 minimization. The parameter p\ represents the characteristic 
time of the flux variation. For the days leading up to the flare (MJD 55291-55295, April 5-9), the 
minimal x 2 /dof value of 0.3/2 is obtained for p\ = 2.87 days. The error bars of the fit parameters 
Po and pi are determined by finding the parameter ranges with x 2 between x m in an d Xmm + 2.30, 
where Xmin ^ s ^ ne smallest x 2 value. The same x 2 calculation is repeated for data from the peak flux 
onward. The details are presented in Table [TJ For the period between MJD 55296 and 55304 (April 
10-18), the exponential decay time is 1.12+0 26 days. An even shorter decay time of 0.90lo'i5 days 
is obtained by restricting the fit to the period between MJD 55296 and 55298 (April 10-12). To 
investigate the possibility of a second flare between MJD 55299 and 55301, a constant-flux fit is 
applied to data points between MJD 55298 and 55304 (April 12-18). The xV dof value of the 
constant-flux fit is 9.6/4 with a corresponding x 2 probability of 0.05. In spite of this low confidence 
level for the constant-flux hypothesis, there is nevertheless insufficient evidence to confirm the 
presence of a second, separated flare component around MJD 55299-55301 (April 13-15). 



3.2.2. Spectral Analysis 

Figure U] shows the spectra measured during the rising period between April 5 and 8 (MJD 
55291-55294), during the peak on April 9 and 10 (MJD 55295-55296), and during the falling period 
between April 11 and 15 (MJD 55297-55301). Power-law fits of the form $ = $ (E /TeV)~ r are 
applied to all three periods, and the corresponding power-law fit parameters are listed in Table EJ 
The spectral index of the peak period differs from that of the falling period by 2.2 a, and from that 
of the rising period by 1.3 a. The peak period has the hardest spectrum of all three periods. 

A hardness ratio (HR) test is also applied to investigate further the possibility of spectral 
variability between these three different periods. The HR may provide more sensitivity as it is 

Table 1. x 2 minimized parameters of the April flare light curve in Figure [3] (Fit function 
$ = po eC" 55290 )/^ 1 ). The error bars of po and p\ are statistical only. 



Period (MJD) 


X 2 /dof 


X 2 Probability 


Po (cm 2 s 1 ) 


Pi (days) 


55291-55295 
55296-55304 
55296-55298 


0.3/2 
23.7/6 
2.1/1 


0.88 
6.0 x 10~ 4 
0.15 


x 10~ 12 
3.61™ 4 x 10- 9 
lA8 + iH x 10- 8 


9 C7+1.65 
z -°' -0.99 

-(0.90^; 22 ) 
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Fig. 3. — Fits to the 2010 April VHE gamma-ray light curve of M87 leading up to the flare and 
trailing the flare, with fit errors included and shown as shaded regions. The exponential timescale 
is 2.87lg 99 days f° r t ne rising flux portion, and 0.901q 15 days for MJD 55296-55298 segment of the 
falling flux. 



Table 2. Spectral power-law fit parameters and hardness ratios for the three periods of the M87 
flare in 2010 April: Rising, Peak, and Falling. Errors given are statistical only. 



Periods 


MJD Date 


Flux Normalization 


Spectral Index 


X 2 /dof 


Hardness Ratio 






Constant <I>o 


r 










10~ 12 (cm^s^TeV -1 ) 








Rising 


55291-55294 


1.92 ±0.42 


2.60 ±0.31 


4.1/4 


0.35 ±0.12 


Peak 


55295-55296 


4.71 ± 0.29 


2.19 ±0.07 


4.3/5 


0.28 ±0.03 


Falling 


55297-55301 


1.10 ±0.16 


2.62 ±0.18 


5.2/4 


0.10 ±0.04 
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Fig. 4. — Spectral measurements during three periods: leading up to the flare (MJD 55291-55294, 
April 5-8), peak of the flare (MJD 55295, 55296, April 9 and 10), and trailing the flare (MJD 
55297-55301, April 11-15). The lines are power-law fits to the data, with the values for the flux 
normalization constant and the spectral index given in Table [2j 
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obtained in a straightforward way from the energy distribution of the excess events, whereas the 
spectral index calculation requires multiple binning and fitting of the data. The HR used here is 
defined as the ratio of the integral flux in the energy range 1-10 TeV to that in the range 0.35- 
1 TeV. For the rising period, HR = 0.35 ± 0.12; for the peak period, HR = 0.28 ± 0.03; and for the 
falling period, HR = 0.10 ± 0.04 (see Table [2|). The HR for the peak period is found to be larger 
than that for the falling period with a statistical significance of 3.6 a, compared to the 2.2 a for 
the spectral index difference of the same time periods. The increased significance may result from 
a higher sensitivity of the HR to spectral variability. However, we also note that the HR for the 
rising period is 2.0 a larger than that for the falling period while the spectral indices from these 
periods are identical. This may be a result of statistical fluctuation due to the poor statistics of 
the rising period spectral measurements. 

Figure [5] shows the spectral index (r) plotted against flux normalization constant ($o) for the 
2010 Ap ril flare spectra (open circles), together w i th archival VHE gamma-ray spectra f rom 2004 



onward (Aharonian et al 



2006 



Acciari et al. 



2008 



Albert et al 



2008 



Acciari et alJl200fl . 12010 ). A 



constant-flux fit to the 2010 April flare flux-index data yields a x probability of 0.05. A linear fit of 
the form T = po +pilog 10 $o yields a \ 2 probability of 0.67, with the parameter p\ = —0.72 ± 0.30, 
which is 2.4 a away from zero. Although the fit may suggest a possible correlation between the 
spectral index and the flux normalization constant, the data do not provide definitive evidence for 
spectral variability during this flaring episode. Using all the flux-index data available since 2004, 
a constant-flux fit yields a x 2 probability of 0.26, while a linear fit yields a x 2 probability of 0.52 
with pi = -0.27 ±0.13. 



4. Discussion 



VERITAS first detected M 87 in 200 7 in a low state, with emission at ~ 2% of the Crab Nebula 
flux above 250 GeV (j Acciari et alJhoOsh . In 2008, VERITAS detected flaring activity up to 10% 
of the Crab Nebula flux above 250 GeV durin g a joint monitoring campaign in which correlations 
between VHE gamma rays, X-ray, and radio (jAcciari et al.1 12009 ) were found. In 2 009, M87 was 
obser ved to be in a low state again at ~ 1% of the Crab Nebula flux above 250 GeV (jAcciari et al 



2010 ). In 2010 April, VERITAS observed the brightest emission ever seen from M87, with a flux 
up to 20% of the Crab Nebula flux above 350 GeV . In comparison to previous constraints from 



past flares (jAharonian et al.ll2006l : I Albert et al 



2008 



Acciari et alJl201(J ). the 2010 VERITAS data 



set yields the fastest exponential flux-changing time (0.90l ^ 5 days) ever observed for M87. This 
time constraint gives a new upper limit on the emission region size that is lower than those derived 
from previous observations. Using the exponential decay time, the emission region size has radius 
R < Rvar — 5 c At = 2.3 x 10 15 5 cm ~ 1.35 R s , where R s is the Schwarzschild ra dius of the M87 
black hole (= 2GM BH /c 2 w 1.8 x 10 15 cm, with M BH = 6.2 x 10 9 M Q scaled from bebhardt et al 



(|201ll ) to the distance used in this paper) and 5 is the relativistic Doppler factor. As in earlier 
findings, this may point to the black hole vicinity as the actual origin of the VHE radiation. While 
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Fig. 5. — Spectral index vs. flux normalization constant using spectra from the three periods 
(rising, peak, and falling) and archival spectra from 2004 onward. The dashed line represents a 
constant fit with a x 2 probability of 0.26, and the dash-dotted line represents a linear fit of the 
form r = pq + pilog^^o with a x 2 probability of 0.52. The values of po and p\ are given in the 
text. 
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an increased X-ray flux from the nucl eus seems to support thi s hypothesis, no increase of the radio 
flux from the nucleus could be found ([Abramowski et alJl201lT). in contr ast to the contemporaneous 
radio and VHE gamma-ray flares observed in 2008 (jAcciari et al.l 12009 ). 



Another notable characteristic of the 2010 flare is the large difference between the rise time 
and the decay time of the flux, a feature which has not been seen in previous flares. Since previous 
VHE flares in 2005 and 2008 were not sampled at a comparable accuracy and their onsets were not 
as well defined as the 2010 flare, this is the first M87 VHE flare that allows the determination of 
the rise and fall times. The shape of the 2010 flare also seems less erratic as compared to the earlier 
flares, which could point to a different production mechanism. However, given the lower statistics 
of the earlier flares, this is difficult to quantify and requires future observations to disentangle. 



From a compilation of multiwavelength data sets spanning decades, I Wagner et al.l (120091 ) pre- 
sented a spectral energy distribution (SED) of M87, a long with hadronic and leptonic models. 
The hadronic synchrotron-proton blazar (SPB) model (jReimer et al.l |2004j) suggests gamma-ray 
emission from synchrotron radiation by protons or by muons and pions. However, the SPB model 
SED produced using archival data before 2004 shows a steep drop-off at TeV energies that is not 
compatible wit h the spectra obtain ed from the 2010 data set or with any previous VHE spectral 



measurements. IBarkov et al.l (|2010l ) proposed a scenario where a red giant star, with an envelope 
of loosely bound material, interacts with the base of the jet. VHE gamma rays are produced near 
the supermassive black hole via proton-proton interactions between the jet and the red giant cloud. 
The gamma-ray light curve produced from this model shows an exponential increase/decay time 
of ~ 1 day, identical to the decay timescale obtained from the 2010 data. However, the model gave 
no prediction on the VHE spectrum for comparison with our data. 

There are also several leptonic jet models with different geomet ric structures that can explai n 
the VHE gamma-ray emission, such as the decelerating j et model bvlGeorganopoulos et al.l (|2005l ). 



the multi -blob synchrotron-self-Compt on (SSC) model by lLenain et al.l (|2008l ). and the spine-sheath 



model by lTavecchio fc Ghisellinil (J2008J). The SED solutions obtained from these models can explain 
the observed rapid variability and match the VHE spectra well in both low and flaring states from 
2005 and earlier. Looking into some of these models in more detail, the model para meters can be 



adjust ed to account for recent measurements. In the case of the multi-blob scenario, iLenain et al 



(|2008l ) showed that for the case of M 87, the model spectrum hardens with decreasing magnetic field. 
In order to keep the size of the VHE emitting region of the order of the Schwarzschild radius, the 
local value of the magnetic field should be < 0.01 G. The spine-sheath model, however, seems to face 
difficulties in achieving a harder spectrum due to absorption of TeV photons from interacti ons with 
the optical-IR photons from the spine. As pointed out by iTavecchio k, Ghisellinil (120081 ). severe 
gamma-ray photon absorption can be alleviated by increasing the emission region size, which would 
decrease th e absorption optic al depth. However, this would be limited by the observed short-term 
variability. Abdo et al. (j2009l ) fitted a homogeneous one-zone SSC model using 2009 VLBA radio, 
Chandra X-ray, and Fermi-LAT measurements when M 87 was observed to be in a low state from 
radio to VHE gamma rays. A contemporaneous spectral measurement in the VHE range was not 
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possible due to low statistical significance (jAcciari et al.ll20ld ). but compared to archival low-state 
VHE measurements, the one-z one SSC model seems to und erest imate the VHE gam ma-ray flux by 
almost an order of magnitude. Georganopoulos et al. ( 2005 ) and Lenain et al. ( 20081 ) demonstrated 
that one-zone homogeneous models are unlikely to reproduce the observed VHE spectrum. 



Giannios et al.l (|2010l ) presented a scenario where minijets are formed within the jet due to flow 



instabilities. These minijets move relativistically with respect to the main jet flow. VHE gamma 
rays are produced from the interactions between the minijets and the jet, and are beamed with 
large Doppler factor when the minijets are aligned with our line of sight. The minijets model SED 
is compatible with the 2010 data. A satisfactory solution for the hi gh state observed by VERI- 
TAS in 2010 is also possible within the magnetosphere model (e.g.. iNeronov fc Aharonianl 120071 ; 



Rieger Aharonian 



2008 



Vincent fc Lebohec 



ctosp 

2010I : 



Levinson Riegei 



2011 



The magnetosphere 

model is dependent on the injected plasma which suggests that a vacuum gap with a large electric 
field that is capable of accelerating electrons to very high energies may be formed during a period 
of low accretion rate. 

We cannot discriminate between different leptonic models based on this VHE data alone. The 
spectral change with flux level would serve as an important input for the modeling once it is 
confirmed by a second flare. Leptonic models tend to predict a more direct correlation between 
X-ray and VHE gamma ra ys. For the 2010 flar e of M87, extensive follow-up observations of 
the VHE gamma-ray flare (lOng Sz Mariottil l2010l ) were carried out in X-ray, optical, and radio 
wavebands. The result is a much more complete sampling across different energy bands than in 
the case of previous M87 flares, providing a data set that will help to con strain the emission region 
and the radiative processes involved. A separate, upcoming publication (jAbramowski et al.ll201ll ) 
will present the multiwavelength result, which spans 16 decades of energy. 
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